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Abstract 

We outline the basic ideas involved in a recently proposed [17] derivation of a gauge 
theory for underdoped cuprates in the "spin-gap phase" , performed essentially step 
by step starting from the t — J model, considered as a model Hamiltonian for the 
Cu02 layers. The basic tool is the £7(1) x SU{2) Chern-Simons bosonization, to 
which it is dedicated a somewhat detailed discussion. The basic output is a "spin- 
gap" not vanishing in any direction and an antiferromagnetic correlation length 
proportional to the inverse square root of doping concentration, in agreement with 
data deduced from the neutron experiments. The model also exhibits a small half- 
pocket Fermi surface around (±7r/2, ±7r/2) and a linear in temperature dependence 
of in-plane resistivity in certain temperature range. 
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1 . High T c Cuprates and the t-J model 

A common structural feature of high T c cuprates is the presence of electronically 
active CuOi layers, alternating with (insulating) block layers along the crystalline 
c-axis. In the CuO^ layers of undoped materials the 3d shell of copper has a 
hole (primarily in the highest energy 3d x i_ y i orbital) while the 2p shell of the 
oxygen is filled. The spin ^ moments of the Cu are antiferromagnetically ordered 
at low temperatures and a strong on-site Coulomb repulsion acts in the 3d orbital, 
inhibiting double occupation. As the materials are doped, holes (or electrons) are 
introduced in the Cu02 layers. 

In terms of doping concentration (8) and temperature (T), a "typical" phase di- 
agram is drawn in fig.l (patterned on La2-sSrsCuO^ compounds). It exhibits 
an antiferromagnetic (AF) insulating phase near 5 = (for sufficiently low tem- 
peratures), a superconducting (SC) phase for an intermediate doping (e.g. 0.1 
< 5 < 0.2 for LaSrCuO compounds). The materials with doping concentration 
exhibiting highest T c are called optimally doped; cuprates with lower or higher dop- 
ing concentration are called underdoped or overdoped, respectively. At optimal 
doping, in the phase diagram above the SC region there is a region characterized 
by an anomalous metallic behaviour (e.g. linear in T in-plane resistivity p a b [1]; 
anomalous spin lattice relaxation rate ~ tst> a ~ 1 [2]; large 2D Fermi surface 
whose volume is consistent with the Luttinger theorem [3]). Moving towards the 
underdoped region there is a crossover to the so-called "spin-gap phase" , exhibit- 
ing distinctive phenomena (e.g. a minimum of p a b(T) for small enough 6 [4]; a 
maximum in — ^ a t low T [5]; in some materials small half-pocket like 2D Fermi 
surface around (±f, ±^)[6]). In the overdoped region the materials appear to 
show an essentially "normal" metallic behaviour. 

In this paper we will be interested in the underdoped (spin-gap) region. A model 
Hamiltonian for the Cu02 planes (at low doping) has been proposed by Zhang 
and Rice [7], following a suggestion by Anderson [8], roughly on the basis of the 
following considerations (see [9] for a more precise and detailed discussion). The 
holes introduced by doping go primarily into symmetrized O-orbitals around the 
Cu ion and they form a spin singlet with the spin moment of copper (see fig. 2). 
A spin singlet (Cw-hole /O-hole) in one CuO^ has a relevant nearest neighbour 
(n.n.) hopping, since each CuO^ has an O-site in common with the n.n. CuO^. 
The low energy physics is then believed to be dominated by the motion of these 
spin-singlets in the AF background of Cu spin moments. The Hamiltonian pro- 
posed for the system is given by 
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H = P G[J2H -MaCj* + h-C.) + JS % ■ Sj]P G , (1.1) 

<ij> a 

where i runs over the sites of the (square) lattice defined by the position of the Cu 
ions, the sum over a runs over spin indices (spin up=l, spin down=2). In eq. (1.1), 
Ci a denotes a spin ^ fermion operator, Pq is the Gutzwiller projection, eliminating 
double occupation, modelling on-site Coulomb repulsion and the second term is 

— * 

an AF- Heisenberg Hamiltonian, where the spin Si is given by 

Si = c\/-f Ci ,. (1.2) 

The system described by the Hamiltonian (1.1) is called u t — J model". 

According to [9], to match with the physics of high T c cuprates one should take 
J l 

t 3 - 

We analyse the t — J model with a path-integral approach (see e.g. [10]). We 
write the euclidean action, corresponding to the Hamiltonian (1.1), with chemical 
potential fx, at temperature T, in terms of spin ^ (Grassmann) fermionic fields 
Wq,, \&* . The action is given by 



rf3 

J ° ' <ij> 



where the two-body potential Uij is given by 

J +oo i=j 

U ^ = \-{ ij n.n ( L4 ) 

and 6 = fi + J/2. The last two terms in (1.3) correspond to a rewriting of the 
AF Heisenberg term plus a hard-core repulsion replacing Gutzwiller projection. 
[Summation over repeated spin indices is understood and dependence on euclidean 
time x° is not explicitly exhibited; (3 = ^r^f , with ks the Boltzmann constant]. 
For example, the grand canonical partition function of the t — J model is expressed 
in path-integral form as 

E(/3, fj) = J V^V^e~ s ^^ (1.5) 
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2. "Chern-Simons representation" of the t-J model 

A key problem is to find a good Mean Field Approximation (MFA) to analyze the 
t — J model. 

A priori, we can consider many possibilities as starting point for a MFA from a 
general procedure valid in 2D: the "Chern-Simons bosonization" . Let be a 
gauge field of gauge group G; define the (euclidean) Chern-Simons action by 

Sc.sXW) = -±- i J^dx°J dxe^tr [W»d u W p + ^W„W u W p ](x , X*)-, (2.1) 

denote by S(^, ip*, W) the action obtained from S(ip,ip*) by minimally coupling 
ip,ip* to the gauge field W. 

The output of "Chern-Simons bosonization" can be summarised as follows: In 
2D, for suitable choices [11] of the group G and of real coefficient (s) ko, one can 
replace the path-integration over by path-integration over W and new spin 

| fields, x, x*, bosonic or fermionic depending on {G, ko}, substituting the action 
by 

S( X ,x\W) + k G S c . s .(W) (2.2) 

(with suitable b.c. [11]) and the new theory is exactly equivalent to the original 
t — J model. 

Remarks 

1) Although the C.S. bosonization is an exact identity if treated without ap- 
proximations, each one of these "C.S. representations" in terms of x, %*, W can 
be taken as a starting point for a different MFA. As an axample, if we choose 
G = C/(l),fc{7(i) = 1 and x, x* bosonic, one can reproduce [12] the "slave-boson" 
approach [13]. 

2) A similar strategy works for the Fractional Quantum Hall Effect: the system 
in a plateau of Hall conductivity around a filling fraction v = ^rpj,^ integer, 
of the first Landau level appears [14] to have a good MFA in terms of a "C.S. 
bosonization" of the original long-distance action with G = U(l), ku(i) =2£+l. 

Here we sketch how the C.S. bosonization works intuitively, while more details are 
given in the Appendix (for a rather complete discussion see [11]). [In the following, 
we usually omit explicit reference to the conjugate fields, like \E'*]. 
Integrating out the time-component of the gauge field, Wo, appearing linearly in 
the action, one obtains a constraint of the form 
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Jo(x) = ^eoupW^ix), (2.3) 

where W^ v is the field strength associated with (in particular W 12 is the 
"G-magnetic field") and j M is the "G-current" of the matter field. 
As a result a "G-vortex" is attached to every particle described by XiX*i hence 
assigning to them a "G-magnetic charge" . These particles also carry a "G-electric 
charge", since x is minimally coupled to W^. The presence of both "electric" and 
"magnetic" charges implies a Aharonov-Bohm (A-B) effect when these particles 
are exchanged, thus introducing a "phase" for every exchange [15]. If this "A-B 
phase" is trivial (+1), then the statistics of x is unchanged after VF-integration 
and, to match with that of the field x m ust be taken fermionic. If the "A- 
B phase" is -1, then integration over W turns the statistics of a bosonic x m to 
fermionic, thus matching with that of \& (after VF-integration) . 
Furthermore one can prove that the only effect of the C.S. term is to introduce the 
"A-B phase". These arguments give an idea why C.S. bosonization is an exact 
identity. 

To analyze the t — J model we choose x fermionic, 

G = U(l) x SU(2), k u{1) = -2, k su(2) = l. (2.4) 

The gauge field with gauge group U(l) is denoted by 5 M and the gauge field with 
gauge group SU(2) is denoted by V^. Why we made this choice? 
The basic argument is that with this choice, applying a dimensional reduction 
{2D — > ID) and performing a MFA we are able [16] to reproduce the exact results 
on the one-dimensional t — J model in the limit t » J, obtained by Bethe- 
Ansatz and Conformal Field Theory techniques, including the "semionic" (i.e. 
intermediate between bosonic and fermionic [15]) statistics of spin and charge 
excitations. 

3. A formal separation of charge and spin degrees of freedom 

In 2D we separate formally [17] the spin and charge degrees of freedom (d.o.f.) 
of x (the fermion of the t — J model in the chosen "U(l) x SU(2) Chern-Simons 
representation" ) by a polar decomposition: 

Xa = #£a- (3-1) 

In (3.1), H is a spinless fermionic field ("holon") and it is minimally coupled to B 
respecting a local U(l) gauge invariance; E a is a spin \ bosonic field ("spinon"), 
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it is minimally coupled to V respecting a local SU(2) gauge invariance and it 
satisfies the constraint 



Kj^aj = 1 (3.2) 

for each site j. An additional abelian local gauge invariance, called here h/s, 
appears due to the decomposition (3.1) of \. It corresponds to the addition of a 
local phase to H and subtraction of the same phase from E a . 

Remarks 

1) Fields with the same quantum numbers and constraints of those appearing in the 
"slave fermion" approach [18] are obtained from H, S a by performing a Holstein- 
Primakoff transformation [12], but the action derived in the U U(1) x SU(2) C.S. 
representation" differs from the slave fermion one. At this stage it is explicitly 
given [17] by 

rP 

S(H, E, B,V)= / dx°{ V[# *(<9 - iB (J) - S)H 5 + iB (J) 
Jo 

+ {l-H*H j )^ a {d +iV Q {j)) a ^ jP ]+ [{-tH;e*f«*> B H^ ia {Pe^<^> V ) 0/9 E i/5 + 

<ij> 

+/*.c.)+^(l-tf*/^)(l-tf*/^ 

(3-3) 

2) The fermionic field $ itself can be written as a product of two U(l) x 577(2)- 
gauge invariant fields, one constructed out of H and B and one out of E a and V; 
the statistics of both these fields is semionic, as originally suggested by Laughlin 
[19]. 

As usual in the path-integral formalism, to proceed we need a gauge-fixing of the 
local gauge symmetries of the action [20]. 

We gauge fix the U(l) symmetry by imposing a Coulomb gauge on B. Integrating 
out the time component So, appearing linearly in the action, one obtains 

B^B^ + SB^H), ^ = 1,2, (3.4) 

- i f B 

where B^ is a mean field introducing a flux it per plaquette p, i.e. e Ja p = —1, 
and 5B^ is a fluctuation term depending on holon density. 

We gauge-fix the SU (2) symmetry retaining the bipartite structure appearing in 
the ground state at zero doping, i.e. the Neel state of the AF Heisenberg model. 
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For this purpose, we adopt the "Neel gauge" defined by Ej a = a^J (g), \j\ = J1+J2, 
where {31,32) denote the coordinates of the site j. The "spins" E*aE are then 
forced in the Neel configuration, leaving all the SU(2) degrees of freedom in the 
gauge field V, unconstrained. We decompose V into a "Coulomb component" V c 
satisfying 

<9% c = o fx = 1,2 

and gauge fluctuations around it, described by an SU(2) field g. Integrating out 
the time-component Vg, appearing linearly in the action, we obtain an explicit 
expression of V° in terms of g and H, denoted by V^{g, H). 

Following a strategy developed in ID [16], with techniques patterned from a proof 
of the diamagnetic inequality [21], we then find a ^-configuration, g m , optimizing 
the partition function of holons in a fixed g-background [17]. 
For small but non- vanishing doping concentration we find that V{H) = V c {g m , H) 
attaches to the holons a vortex of S'L/ (2)-vorticity ±cr z ^, where a z is the diagonal 
Pauli matrix, and the sign depends on the Neel sublattice where the holon is 
located. Furthermore we have gJVx ' = a^ +1 . We rewrite 

V^(g,H) = V fl (H) + 8V^g,H), // = 1,2, (3.5) 

where the fluctuation term SV depends also on the spin d.o.f. described by g, 
whereas V is independent of them. 

Remark 

At this stage, performing a simple field redefinition, the action of the t — J model 
can be exactly written [17] in terms of H and g as S = Sh + S S: where 

S h = f\x°{J2H*(d - (a^gpogja^U - 5)H 3 
Jo j 

+ E [-tH*e-*U> B+SB H^g t P(e^> * + %^!?)n + h.c.]} 

<ij> 

S s = f dx°{J2(^g i j d g j a^) 11 + 
Jo 

+ E ^(1 - - ^*^) [|(^l^(e^<^> ^ V )^^l) 11 | 2 - i]} (3.6) 

<ij> 
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We now make the first basic approximation neglecting 5B^ and dV^, i.e. the feed- 
back of charge fluctuations on B and of spin fluctuations on V c , still retaining 
the holon dependence of V c . Presumably the main neglected effect is a statistic 
transmutation giving rise to semionic statistics for holons and spinons, as in ID. 
We argue that the statistics here is less relevant than in ID because we expect in 
2D the formation of a bound state with the quantum numbers of the electron, due 
to gauge fluctuations discussed later on. 

4. Low energy effective action for spin degrees of freedom 

To derive a low-energy effective action for the spin d.o.f., we first rewrite g in CP 1 
form, introducing a spin | field b a through 

* (i:: #) <"> 

with the constraint b* a bj a = 1 for every site j. 

Subsequently, we apply to b a the standard treatment of AF systems splitting it 
into an AF component, described by a spin | complex boson field z a , a = 1, 2 and 
a ferromagnetic component which is then integrated out [22]. 
The continuum action is then given in CP 1 form by 



S(z,A)= fdx° J dx{\(d -A )z a \ 2 + v 2 8 \(d^-A^z a \ 2 + V\H)z* a z a }(x x), 

(4.2) 

where v s is the "spin velocity", J-dependent, and A is the standard "Hubbard- 
Stratonovich" gauge field of the CP 1 models. [In the derivation we implicitly 
assumed that the CP 1 model is in the symmetric phase; this is self-consistent 
with the behaviour discussed later on]. 

This action describes spin waves Q = z*az interacting with the vortices appearing 
in V 2 , centered at the holon positions. We evaluate approximately the effect of V 2 
averaged over holon positions at fixed density 5. One can argue that this treatment 
can be justified for J « t, because the holon appears to develop a large effective 
mass, due to the coupling to soft spin fluctuations [23]. 
The output is an average 

(V 2 )~-5\n8. (4.3) 

Substituting in (4.2) V 2 by < V 2 > produces a mass for the spin waves, suggesting 
that the system should exhibit short-range AF, with a correlation length £af ~ 
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(— 51n5)~5(~ S~h for 5 not too small). This result is in agreement with data 
deduced from neutron experiments in La2-sSrsCu04, where a fit for £af was 
suggested [24] in terms of 5~^, as shown in fig. 3. 

Remark 

If S(z, A) would be the full action, the spinons z would be logarithmically confined 
by the Coulomb interaction mediated by A, due to the massive nature of z, but 
coupling with holons yields deconfinement. 

5. Low energy effective action for charge degrees of freedom 

Introducing a flux tv per plaquette, the mean field B induces a partition of the 
lattice in two Neel sublattices (here denoted | and j) and it converts the spinless 
fermion H into 2-components Dirac-like fermions of two-species: 

each component of them being supported on a Neel sublattice, as indicated by the 
subscript. The vertices of the cones of the "Dirac" energy-momentuum dispersion 
relation are centered at the four points (±^,±^-) in the Brillouin zone. (This 
phenomenon is standard in the flux phase; see [22,25]). Up to a short-range term, 
the continuum action for these fermions is given by 



,/3 , 2 
S(^,A)= / dx° / dx{Y^^ r \ lQ {d Q -8-e^A ) 
Jo J r=1 

where t is a renormalized hopping parameter and - 
[We adopt the standard notation 

7o = cr 2 7m = ( ff v> a x) 4 = 7m^m $ = ^*1o }■ 

The components of the fermion fields in (5.2) supported in the two Neel sublattices 
have opposite charge with respect to A, which is the same gauge field appearing 
in (4.2). It can be traced back as the gauge field of h/s gauge invariance. 
With respect to the action discussed in [25] , a crucial difference is the appearance 
of the 7o5 term. Neglecting at first A, this term produces a finite Fermi surface for 
the gapless components (i/j^^i/j^). The other two components have a gap and 



-t(^-e (r) 4)]^ (r) }(x°,x), 

(5.2) 

■(i y.-M-o 1 !)■ 



9 



they mix with the gapless ones, due to the presence of the (non-diagonal) Dirac 
7 matrices. 

This mixing is expected to produce a reduction of the spectral weight in the outer 
part in the reduced Brillouin zone scheme. 

In a similar lattice model (with a twist of statistics between holon and spinons) the 
shape of the Fermi surface for the "electron" deduced in mean fields is half-pocket 
like [26], showing a qualitative agreement with the F.S. deduced from ARPES 
experiments in underdoped cuprates [6] (see fig. 4). 

6. Final comments 

We can summarize the large distance behaviour of the model in U(l) x SU(2) C.S. 
representation (with SB = SV = 0), in terms of a system of spinons, z a , whose 
dynamics is described by a CP 1 model with mass term m ~ 5^ (the main novelty) 
and a system of Fermi liquid holons (obtained by integrating out the gapful Dirac 
modes) with ~ tS, interacting via an abelian gauge field A. 
The U(l) effective action for A then exhibits the basic features, such as the Reizer 
[27] singularity in the propagator of the transverse component of A and the ex- 
istence of a characteristic energy scale for spinons, that permits one to apply the 
ideas of [28] to deduce a linear in T in-plane resistivity for certain temperature 
range. 

A more careful study of in-plane resistivity in the present model is in prgress [29]. 
The analysis presented in this paper applies only to small doping concentration 8, 
since we derived under this condition (see eq. 3.5-4.3) the form of V used in the 
treatment, and crucial to get a mass term for spinons. Hence the above description 
is applicable to underdoped materials, as explicitly stated in Sec.l. The extension 
to other regions of the phase diagram is an open problem for further investigation. 

Appendix: Chern-Simons bosonization 

We start by rewriting the grand canonical partition function S(/?, fj,) of the t — J 
model in terms of the canonical partition function Zjy at fixed number of fermions 



N 

Following a standard treatment of Feynman path-integral in the continuum [10], 
adapted to the lattice, we can express Zjv as a sum over virtual trajectories of 



N: 




(Al) 
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fermions on a lattice with imaginary time. These trajectories start at the imagi- 
nary time x° = at a set of lattice sites {ji, j^} and they end at the imaginary 
time x° = (3 at the same set of sites permuted, {j a (i)i ■■■■> ja{N)}i where a is a 
permutation (see fig. 5). Each term in the sum has a weighting factor {—iy^\ 
where e(cx) is the number of exchanges in a. By omitting this factor one repro- 
duces the canonical partition function, Z b N , of a fictitious "bosonic t — J model" 
obtained by substituting fermions with bosons in the t — J model. The hard-core 
constraint forbids any intersection among the trajectories of the N fermions and, 
by periodicity in imaginary time (the planes x° = and x° = (3 are identified, 
see [10]), every set of virtual trajectories defines a link C (see fig. 5). The factor 
(_l)e(o-) i s a topological invariant associated to that link, i.e. it does not change 
under an arbitrary deformation of the link performed without allowing intersec- 
tions. According to a general theory [30] every topological invariant of links can 
be represented as the expectation value of a "Wilson loop" supported on the link, 
in a gauge theory with Chern-Simons action S c , s ,. 

More precisely, let denote a gauge field of gauge group G and ko a real 
coefficient. Then, for a suitable choice of G and kc (for explicit conditions, see 
[11]), we have 

(_!)*(*) = J VWe -k G S c . s .(W )Tr p{e rJ c W^^ (A2) 

where £ is a link associated with a set of virtual trajectories whose end points are 
obtained by a permutation a from the initial points. In (A.2) the last factor is 
the "Wilson loop", where Tr denotes the normalised trace and P(-) denotes the 
normal ordering, which amounts to the usual time ordering T(-) for a fictitious 
"time" parametrising the link. Furthermore, the normalization is chosen such that 
JVWexp{-k G S c . s .(W)} = l. 

Let us denote by Zj^(W) the canonical partition function of the modified "bosonic 
t — J model" with the boson field minimally coupled to the gauge field W. It 
is easy to verify (see e.g. [30] for the abelian and [12] for the non-abelian case) 
that the dependence of W in Z b N (W) appears in the form of linear combinations 
of Wilson loops on the links associated with virtual trajectories. Therefore the 
identity (A.2) implies that 

Z N = J VWe- kGScs(w) Z b N (W). (A3) 

Equation (A. 3) is a "bosonization formula" and it holds for all the couples of G 
and ko satisfying (A.2). If we plug (A. 3) into (A.l), we obtain 
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E = J VWe- kaSc ° Wz b (W), (AA) 
where E b (W) is the grand-canonical partition corresponding to Z b N {W), given by 

E b (W) = J T> X T>x*e- kGS ( x > x *' w \ (A5) 

with x a bosonic field. Combining equations (A. 4) and (A. 5), one obtains the 
"Chern-Simons bosonization" discussed in the text. It is obvious that introducing 
further Chern-Simons gauge fields we can turn x into a fermion, following a similar 
procedure. 

Remarks 

1) A warning: in the above discussion subtle points like boundary condition and 
"framing" are completely ignored, see [11,30]. 

2) To have a more concrete feeling about formula (A. 2), let us consider the simple 
case G = U(l),kjjm = 1, a choice for which the formula holds, and, due to the 
abelian nature of G, path-ordering is not needed, and the trace Tr is trivial. 

Let us denote by Ti^dx^ a "singular <5— like current" supported on a surface E, 
whose boundary is given by the link C. Then one easily verifies that 

e i f c W„dx» =e iJ W^"d^ p d s x_ ( A gj 

We now compute 

J VWe^ 5 ^ p w^w p d 3 x e i J w^d„v p d 3 x (A.7) 

by a change of variable. Shifting — > + 2-zrE^ and "completing the square" , 
one obtains 



/ 



The integral in the exponent of the last term in (A. 8) receives a contribution only 
when £ crosses the surface E, i.e. for every crossing appearing in £, i.e. for every 
exchange of the particles in the virtual trajectories, so that its value is given by 
(-l) e(fj) , q.e.d. . 
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Fig. 1. Presumed phase diagram. 
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Fig. 2. Zhang-Rice singlet. 




Fig. 3. Instantaneous spin correlation length vs temperature in La2_§SrgCuC>4. 
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Fig. 4. Fermi level crossing from two Bi2212 samples of differing oxygen 
content (up) compared with the momentum distribution of "electrons" found 
in [26] (down). 
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Fig. 5. Heavy lines describe a set of virtual trajectories for N=4; L is the link 
obtained identifying the x®=fi and x®=0 planes. 



16 



